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Abstract 

We determine the fermion mass dependence of Euclidean finite volume partition functions 
for three-dimensional QCD in the e-regime directly from the effective field theory of the 
pseudo-Goldstone modes by using zero-dimensional non-linear er-models. New results are 
given for an arbitrary number of flavours in all three cases of complex, pseudo-real and real 
fermions, extending some previous considerations based on random matrix theory. They 
are used to describe the microscopic spectral correlation functions and smallest eigenvalue 
distributions of the QCD3 Dirac operator, as well as the corresponding massive spectral sum 
rules. 



1 Introduction and Summary 



The e-regime of quantum chromo dynamics in four spacetime dimensions (QCD4) [1, 2] provides 
one of the few examples wherein exact nonperturbative results can be derived from QCD. It 
corresponds to the infrared sector of the Euclidean Dirac operator spectrum which is related to 
the mechanism of chiral symmetry breaking. The e-expansion of chiral perturbation theory in 
this case, which describes the dynamics of the Goldstone modes [3], is taken with respect to the 
small quantity e = where L is the linear size of the given spacetime volume. If L is much 
larger than the QCD scale, then spontaneously broken chiral symmetry dictates the interactions 
and the leading contributions come from the zero-momentum pseudo-Goldstone fields. This 
regime represents a domain wherein a detailed understanding can be achieved of the lowest- 
lying Dirac operator eigenvalues, which dominate the physical low-energy observables, and in 
which numerical computations in lattice QCD are possible [4]. Most of the recent understanding 
of the e-regime of QCD comes from its relation to universal random matrix theory results [5, 6]. 

In this paper we will study the infrared behaviour of the Dirac operator spectrum in three- 
dimensional QCD (QCD3). In this case the analog of chiral symmetry breaking is played by the 
spontaneous breaking of flavour symmetry [7] . The finite volume gauge theory is connected to 
three universality classes in random matrix theory, the unitary ensemble (fundamental fermions 
with N c > 3 colours) [8, 9], the orthogonal ensemble (fundamental fermions with colour group 
SU(2)) [10] and the symplectic ensemble (adjoint fermions with colour group SU(N C )) [11]. The 
random matrix theory representation of the effective QCD3 partition function in the e-regime 
has lead to a number of conjectures regarding the derivation of the spectral properties of the 
three-dimensional Dirac operator from joint eigenvalue probability distributions in the pertinent 
matrix model [8, 12]. The patterns of flavour symmetry breaking in all three instances have been 
identified [8, 10, 11] and supported by Monte Carlo simulations [13], and the microscopic spectral 
correlation functions have been calculated [6],[14]-[18]. 

To derive nonperturbative analytic results for the smallest eigenvalues of the Dirac operator 
it is necessary to have a complete proof that the universal random matrix theory results coincide 
with the low-energy limit of the effective field theory. From the field theoretic point of view [2] 
these observables should be computable entirely within the framework of finite volume partition 
functions which are generating functionals for the order parameter for flavour symmetry break- 
ing [15]. For QCD3 this problem has only been thoroughly investigated for the unitary ensemble 
by using the supersymmetric formulation of partially quenched effective Lagrangians [19], the 
replica limit [20], and the relationship [21] between finite volume partition functions and the 
r- function of an underlying integrable KP hierarchy [22, 23]. The goal of this paper is to derive 
finite volume partition functions in all three cases mentioned above and to use them to explore 
spectral characteristics of the QCD3 Dirac operator directly in the low-energy effective field 
theory. 

In all three ensembles the effective field theory is described by a non-linear u-model in a 
maximally symmetric compact space [8, 10, 11]. These target spaces fall into the classical 
Cartan classification of symmetric spaces, and confirm part of the Zirnbauer classification of 
random matrix universality classes in terms of Riemannian symmetric superspaces [24]. This 
enables us to exploit the well-known geometrical properties of these spaces [25, 26] (see [27] 
for a recent review in the context of random matrix theory) and hence explore the differences 
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between the three classes of fermions in QCD3. For example, we will find that the finite volume 
partition functions for adjoint fermions are drastically different in functional form from those for 
complex and pseudo-real fermions, and this difference can be understood from the form of the 
Riemannian metrics on the corresponding symmetric spaces. A similar geometric feature also 
exhibits the difference between even and odd numbers of quark flavour components. In the latter 
cases the theory is plagued by a severe sign problem which is rather delicate from the effective 
field theory point of view [8]. We derive new results in all three instances for the finite volume 
partition functions with odd numbers of flavours. In the case of the unitary ensemble, we use 
them to derive a new simplified expression for the microscopic spectral correlators, extending 
the known effective field theory results for the even-flavoured case [19]. We also study some 
of the spectral sum rules which constrain the eigenvalues of the Dirac operator [8, 2, 28], and 
briefly describe its smallest eigenvalue distribution functions [29]. 

On a technical level, we present a complete and thorough analysis only for the case of complex 
fermions. In the cases of pseodo-real and real fermions, we complete the derivations only in the 
cases of completely degenerate mass eigenvalues, as there are some algebraic obstructions which 
prevent us from obtaining closed analytical expressions for the finite volume partition functions 
for arbitrary configurations of quark masses. Nevertheless, we present a very general algebraic 
description of the differences and similarities between the three universality classes, and also of 
the distinction between the even and odd flavoured cases, based on a remnant discrete subgroup 
of the flavour symmetry acting on the mass parameters of the gauge theory. This residual 
symmetry is described in general in Section 2, where we also review various aspects of the 
non-linear cr-models and describe the geometry of their symmetric target spaces in a general 
setting that applies also to other symmetric space a- models such as those encountered in QCD4. 
Sections 3-5 then use this general formalism to analyse in detail the three classes of fermions in 
turn. In addition to obtaining new results in the odd-flavoured cases, the main new predictions 
of the effective field theory approach concern the (equal mass) partition functions for adjoint 
fermions, which have a far more intricate functional structure than the other ensembles and 
which involve functions usually encountered in QCD4 [2, 30], not in QCD3. 



2 Finite Volume Partition Functions 

The Euclidean Dirac operator governing the minimal coupling of quarks in the QCD3 action is 
given by 

i# = i<7" (d^ + iA^) , (2.1) 
where cr M , ji = 1, 2, 3 are the Pauli spin matrices 

and are gauge fields of the local colour group SU(N C ), N c > 2. When multi-colour fermion 
fields are coupled to this operator, the three-dimensional field theory also possesses a global 
flavour symmetry described by rotations of the quark field components via elements of some 
Lie group G. If one introduces quark masses in equal and opposite pairs, then the fermion 
determinant is positive definite and the quantum field theory is anomaly free [8]. The addition 
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of such mass terms preserves three-dimensional spacetime parity [16] but breaks the flavour 
symmetry, and thereby permits one to analyse the possibility of spontaneous flavour symmetry 
breaking in QCD3. In this section we will construct the general low-energy effective field theory of 
the associated pseudo-Goldstone bosons in the ergodic e-regime wherein exact, nonperturbative 
information about the spectrum of the QCD3 Dirac operator (2.1) can be obtained. 

2.1 Non-Linear er-Models of QCD 3 

We are interested in the mechanism of flavour symmetry breaking in three-dimensional QCD, 
which we consider defined at fixed ultraviolet cutoff A, while the quantum field theory in the 
infrared is regulated by a three- volume V. The effective field theory for the low- momentum 
modes of the Nambu-Goldstone bosons can be derived by using the three-dimensional version of 
chiral perturbation theory [3]. It is consistent with the flavour symmetry of the original quantum 
field theory represented through transformations of the quark fields by the Lie group G. In the 
ergodic regime, the effective Lagrangian in Euclidean space is given by 



where f n is the decay constant of the pseudo-Goldstone modes, S is the infinite-volume quark- 
antiquark condensate which is the order parameter for flavour symmetry breaking, and M is 
the Nf x Nf mass matrix induced by integrating out the Nf > 1 flavours of fermion fields in 
the QCD3 partition function. When the number of fermion flavours Nf = 2n is even, so that n 
is the number of four-component spinors, we assume that the eigenvalues of M occur in parity- 
conjugate pairs (mi, —mi), mi > 0, i = 1, . . . , n. When Nf = 2n + 1 is odd, there is a single 
unpaired quark of mass tjiq > which breaks parity symmetry explicitly when tuq > and 
radiatively when m = [31]. The higher order terms in (2.3) represent contributions involving 
gluons and confined quarks. The symbol Tr , in the representations that we shall use in this 
paper, will always correspond to an ordinary Nf x Nf matrix trace. The Nambu-Goldstone 
field IA lives in the appropriate vacuum manifold for the given symmetry breaking. We will use 
the Dyson index (3 = 1, 2, 4, which is borrowed from random matrix theory terminology [32]. It 
labels the anti-unitary symmetries of the QCD3 Dirac operator which in turn depends on the 
types of fermions that are present in the original field theory [9]. 

The Goldstone manifold is a coset space 



where H-p is the stability subgroup of the original flavour symmetry group G which leaves fixed 
an Nf x Nf parity matrix V representing an involutive automorphism of G that reflects the 
parity symmetry of the original field theory. The fields in (2.3) may then be parametrized as 



with U G Gp(Nf). In (2.5) the superscript t denotes the appropriate conjugation involution on 
the Lie group G, given by 17* = f7 f for = 2 and = U T for (3 = 1, 4. 

Given the non-linear cr-model (2.3), we would now like to study the region where the zero- 
momentum mode of U dominates. This is the e-regime in which the linear dimension of the 
system is much smaller than the Compton wavelength of the Nambu-Goldstone bosons. The 




(2.3) 



G p (N f ) = G/H r , 



(2.4) 



U = UTU t 



(2.5) 
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spacetime integration over the effective Lagrangian (2.3) then produces an overall volume fac- 
tor V, and the partition function of the non-linear cr-model simplifies to a finite-dimensional coset 
integral. The region of interest is obtained by simultaneously taking the thermodynamic limit 
V — ► oo and the chiral limit M. — > Onj of massless quarks, with the rescaled masses M = V X M 
finite. The effective field theory of the Goldstone bosons is then the correct description of the 
dynamics below the scale set by the mass gap. Because of the tuning of the quark masses to 
zero, the higher-order corrections to this effective theory are exponentially small in the mass 
gap and come from integrating out the heavy states. 

The quantum field theory is thus adequately described by the finite volume partition func- 
tion [8, 10, 11] 



zf f) {M) = j du p {U) exp Re Tr (MUTU 1 ) 

+ \ (l-(-l) qN f) {-1) N P exp ReTr (-MUTU*) 



(2.6) 



where qjy f is the even/odd congruence class of the flavour number Nf given by c/2n = 0, q2n+i = 1> 
and Np is the number of Dirac operator eigenvalues (defined, for example, in a lattice regular- 
ization or in random matrix theory). The second line of (2.6) only contributes when there is an 
odd number of fermion flavours. It arises from the fact that, generally, the QCD3 partition func- 
tion changes by a phase ( — l) N f N P under a parity transformation of the quark masses. When 
Nj = 2n + 1 is odd, the fermion determinant is not positive definite and its sign determines 
two gauge inequivalent vacuum states which both have to be included in the low-energy limit 
of the QCD3 partition function in the e-regime [8]. We will always ignore irrelevant numerical 
constants (independent of M) in the evaluation of (2.6), as they will not contribute to any 
physical quantities and can be simply absorbed into the normalization of the partition function. 
The measure dvp(U) is inherited from the invariant Haar measure for integration over the Lie 
group G. It will be constructed in the next subsection. 

The low-energy effective partition function of the zero-dimensional cr-model (2.6) contains 
sufficient information to completely constrain the low-energy Dirac operator spectrum. For 
instance, by matching exact results from random matrix theory with the low-energy effective 
field theory (2.6), it is possible to express the microscopic limit of the spectral fc-point correlation 
functions 

^(Ai, • • • , \ k ;M) = / J] Tr <5(A ; -i|) + iM)\ (2.7) 

\ i=i I a 

of the QCD3 Dirac operator (2.1) in terms of a ratio of finite volume partition functions, namely 
(2.6) and one which involves (3k additional species of valence quarks of imaginary masses i A; [15, 
22]. In (2.7) the average is taken over SU(N C ) gauge field configurations weighted by the three- 
dimensional Yang-Mills action. Both (2.6) and the microscopic limit of (2.7) depend only 
on the eigenvalues ±//j,i = l,...,n (and no when Nf is odd) of the scaled mass matrix M. 

For (3 = 2, 4 one explicitly has 
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^%,...,a;M) = cf]](iC 1 + M r/]](c i 2 + ^) TIK-t 2 / 

1=1 1=1 m<p 

x ( } 



where Q = V £ A; are the unfolded Dirac operator eigenvalues, and the proportionality constant 

(k) 

Cp is fixed by the matching condition between the confluent limit of the microscopic correlators 
and the macroscopic correlation functions at the spectral origin as 

k 

For (3 = 1,2,4, the finite volume partition functions may be used to compute the normalized 
hole probability E^ f \(;M), i.e. the probability that the interval [— £, (] C K. is free of Dirac 
operator eigenvalues. One has the formula [29] 

£ ™ K . M )-e-JC>f£!V™) (210) 
(<1 ' Z f'»(M) ' ' ' 

from which the smallest eigenvalue distribution is given by 
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Jf"(C;M) = g ^ ' . (2.11) 



Furthermore, by matching (2.6) with the formal expansion of the full QCD3 partition function 
Z^ l \M) = lf{&et(lp 2 + miy det(^ + im o r/ p M (2.12) 

with pi = p2 = 1 and ^4 = 5, one can express the mass-dependent susceptibilities in terms 
of derivatives of the finite volume partition functions with respect to one or more quark mass 
eigenvalues [28]. This leads to a set of massive spectral sum rules once (2.6) is known explicitly. 
The simplest such sum rule expresses the mass-dependent condensate for the flavour symmetry 
breaking in the case of equal quark masses as 

/y 1 \ _ 1 01nZ^V,...,„) 

In the massless case p — ► 0, the left-hand side of (2.13) can be computed directly from the coset 
integral (2.6) [8]— [11] and is given as a function of the Dyson index (3 by 

\(>o^ I p 2{N f -\) y d + 1 
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for even values of Np. The formula (2.14) generalizes the one given in [10, 11, 17, 33] to include 
the case of odd Nf. For f3 = 2 it was first obtained in [8]. Note that for an odd number of Dirac 
operator eigenvalues, the finite volume partition function (2.6) is not positive and vanishes when 
the unpaired fermion is massless, i.e. /j,q = [16]. It is therefore not suitable to describe physical 
quantities such as chiral spectral sum rules, while it can appear radiatively [15]. 

2.2 Parametrization of the Goldstone Manifold 

We will now describe how to evaluate the non-linear a-model partition functions (2.6). For this, 
we shall first review some elementary Lie theory that will be of use to us and which will fix some 
notation. We assume that the quark flavour symmetry is represented by a compact, connected, 
simple Lie group G, which we take to be realized by Nf x Nf invertible matrices U. Let t be a 
Cartan subalgebra of the Lie algebra of G, so that T = exp(it) is a maximal torus of G. Any 
element X £ Q can be taken to lie in the Cartan subalgebra t without loss of generality, since 
this can always be achieved via rotations by elements of G if necessary. 

After performing such a rotation in the given Cartan decomposition of G, the Jacobian for 
the change of Haar integration measure is given by a Weyl determinant Aq{X) of G which is 
determined through the product 

A G (X) = H(a,X) , (2.15) 

where (a, X) are the positive roots of G evaluated on the Cartan element X. Its most useful 
form is obtained by choosing an orthonormal basis e*j, i = l,...,dimT of vectors in weight 
space, e*j • e*j = Sij. Then the root vectors a can be identified with dual elements 

dimT 

a v = > ( 2 - 16 ) 

i=i 

while 

dimT 

X= ^x,^ (2.17) 
i=i 

where Hi are the generators of t in the matrix basis corresponding to the orthonormal weights. 
By identifying (2.17) with the weight vector x v = Yl,i x i^ii the pairing in (2.15) may then be 
written as the inner product 

dimT 

{a, X) = a v • x v = ^2 en Xi . (2.18) 
i=i 

With U = exp(iX), the suitably normalized invariant Haar measure on the Lie group G is 
given by 

dimT 

[dU] = Yl n sinh2 ( a > ln[/ ) • ( 2 - 19 ) 
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G 


Positive Weights 




Dynkin Diagram 


Weyl Group 


U(r) 

3=2 


1 < i < j < r 


1 < i < r 


o— o— o • • • o-o 




Sp(2r) 
3=1 


e*i ± e j , 1 < i < j < r 
2e*j , 1 < i < r 


1 < i < r 


CE3— •• • • • • 


S r X (Z 2 ) r 


S0(2r + 1) 

3=4 


Pa ~h P a 1 <C i <C 1 <C r 

e*i , 1 < i < r 


-&i-\-r,i-\-r •> 

l<i<r 


a — n o • • • n n 


Sr x (ZoT 


SO(2r) 

8=4 


6-i i 6j , 
1 < « < j < J" 


1 < i < r 


o—o ■ ■ ■ o-^Lo 


S r X (Za^" 1 



Table 1: T/ie flavour symmetry groups G of rank r = dimT > 2 and their corresponding Dyson 
index (3. The vectors e*j, i = 1, . . . ,r form an orthonormal basis of weight space, while E^ u is 
a corresponding orthonormal basis of matrix units (E^^xp = <^a &vp- The Dynkin diagrams all 
have r nodes. The symmetric group S r permutes the eigenvalues Xi of a given Cartan element 
X £ i, while the cyclic group 7L*2 reflects them. 

Its infinitesimal form induces a measure on the Lie algebra q as 

dimT 

[dX] = [] dxi A G (X) 2 . (2.20) 
i=i 

For the Lie groups that we shall deal with in this paper, the Cartan subalgebra t is represented 
in the basis (2.17) by Nf x Nf diagonal matrices. 

The residual flavour symmetry after such a transformation is given by the discrete Weyl 
group Wq of G which consists of the inequivalent transformations on t — > t given by the adjoint 
actions 

X\ — >X w = wXw~ 1 , weG. (2.21) 

Geometrically, Wq is the group of automorphisms of the root lattice of G corresponding to Weyl 
reflections. To each w 6 Wq we may associate a sign factor sgn(u>) = ± 1 which is the even/odd 
Z2-parity of the Weyl element according to its action (2.21) on the diagonal elements of the 
Cartan matrices. For the partition functions that will be of interest to us in this paper, the 
relevant group theoretical data are summarized in Table 1. 

We are now ready to construct the measure dz^g(i7) for integration over the Goldstone man- 
ifold (2.4). For this, we shall introduce polar coordinates on the coset space Qp(Nf) [25], which 
in each case will be an irreducible, maximally symmetric compact space falling into the Cartan 
classification. Let 

= f)r©f)f (2-22) 

be a Cartan decomposition of the simple algebra q with respect to the stability subgroup Hy = 
exp(if) r ), with f)p = f)p the orthogonal complement of f) r in q with respect to the Cartan- 
Killing form. The condition for (2.4) to be a symmetric space is then 

for , f) r ] C f, r , [l) r , tf] c f)f , , f)f] C f, r . (2.23) 
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Let tr be a maximal abelian subalgebra in the subspace f)p, and let A/"t r be the normalizer 
subgroup of tr in Hr with Ct r C A/t r its centralizer subgroup. As Hr is a compact symmetric 
subgroup of the flavour symmetry group G, every coset element U £ Gp{Nf) can be written by 
means of the adjoint representation of the stationary subgroup Hr as [25] 

U = VRV~ 1 , (2.24) 

where V £ Hr/Ct r while R £ Tr = exp(itr) up to transformations by the adjoint actions 
of elements of the factor group jV^,/Ct r , which can be identified with the Weyl group of the 
restricted root system on the symmetric space Gp(Nf). This new root lattice differs from that 
inherited from the original group G if the Cartan subalgebra t is a subset of f) r , and it is defined 
by completing the maximal abelian subalgebra t r C f)p with the generators in f) r that commute 
with t r to give a new representation of t that lies partly in f)p. The matrix R in (2.24) is the 
radial coordinate of the point U £ Qp(Nf), while V is its angular coordinate. This defines a 
foliation of the Goldstone manifold (2.4) by conjugacy classes under the adjoint action of the 
stability group Hr of the symmetric space. 

The decomposition (2.24) means that every matrix U £ Gp{Nf) can be diagonalized by a 
similarity transformation in the subgroup Hr, and the radial coordinates {r{\ are exactly the 
set of eigenvalues of U. Using the Haar measure (2.19) on the Lie group G, the corresponding 
Jacobian for the change of integration measure by (2.24) can be computed by standard techniques 
and one finds [26, 34] 

dim T r 

dvp(U) = [dV] Yl dr i II |smh(a r ,lni?) , (2.25) 

i=l o r >0 

where lni? = i J2i r i Hi with ri £ [0, ^] for i = 1, . . . , dimTp. The second product in (2.25) goes 
over positive roots of the restricted root lattice on Gp(Nf), and m ar is the multiplicity of the root 
a r in the given Cartan decomposition, i.e. the dimension of the subspace of raising operators 
corresponding to the root a r in the algebra h r [26]. One can further exploit the ambiguity in 
the choice of radial coordinates R labelled by the restricted Weyl group A/t r /Ct r to alternatively 
parametrize the radial decomposition (2.24) in a way more tailored to the evaluation of (2.6). 
The particular details of this parametrization will depend in general on the random matrix 
ensemble in which we are working. As we will see, the radial flavour symmetry group A/t r /Ct r 
of the Goldstone manifold naturally distinguishes the cases j3 = 1, 2 from j3 = 4, and in all cases 
provides a nice algebraic distinction between the cases of an even or odd number of flavours Nj. 

2.3 The Harish- Chandra Formula 

The computation of (2.6) can be compared with the Harish-Chandra formula which computes 
the group integral [35] 

J [dU] exp Tr {XUYIJ- 1 ) = ^^^(y) E s ^ ex P Tr ( XY ^ ( 2 ' 26 ) 

for X, Y £ t. It expresses the group integral as a finite sum over terms, one for each element 
of the Weyl group, and it is intimately related to the Fourier transform of the Weyl character 
formula for the flavour symmetry group G. Its main characteristic is that the integration in 



9 



(2.26) may be rewritten as an integral over a coadjoint orbit, which is a symmetric space to 
which the Duistermaat-Heckman localization formula may be applied [36] . The Weyl group can 
be described as the factor group Wq = Ng/T, where Mg is the normalizer subgroup of fixed 
points of the left action of T on the orbit G/T. In other words, the right-hand side of (2.26) 
is the stationary phase approximation to the group integral, and the Harish-Chandra formula 
states that in this case the approximation is exact. For this to be true it is important that the 
matrices X and Y live in the Lie algebra q of G. 



3 Fundamental Fermions 

The simplest instance is when the quarks are described by complex fundamental Dirac fermion 
fields with N c > 3 colours. Although these models have been the most extensively studied ones 
thus far, it will be instructive for later comparisons and simply for completeness to rederive 
the finite volume partition functions in this case using the formalism of the previous section. 
In matrix model language it corresponds to the unitary ensemble (3 = 2. The global flavour 
symmetry group is then U(Nf), while the mass and parity matrices are given in a suitable basis 
of the vector space C N f by 



M = < 



Mi 



a 



V 



/in - 



/II 



a 



N f = 2n, 



, N f = 2n + 1 



Hn / / 



a 



r 



i 



(T S <g> l n 



Nf = 2n , 
N f = 2n + 1 . 



(3.1) 



The mass matrix explicitly breaks the flavour symmetry, and its pairing into quark masses of 
opposite sign preserves three-dimensional spacetime parity. The Goldstone manifolds are [8] 



Q 2 (2n) = U(2n)/U(n) x U(n) , 
Q 2 {2n + l) = U(2n + 1) / U(n + 1) x U(n) 

of real dimensions dim^(2n) = 2n 2 and dim^2(2n + 1) = 2n(n + 1). 
3.1 Partition Functions 



(3.2) 
(3.3) 



N f = l 

The case of a single massive quark corresponds to formally setting n = in (3.3). In this case 
there is no flavour symmetry breaking, the Goldstone manifold is a point, and there are no 
massless propagating degrees of freedom. The finite volume partition function is then trivially 
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obtained from (2.6) in the form 



cosh/U 
sinh fx 



Np even 
Njp odd 



(3.4) 



Nf = 2n 

The restricted root lattice of the symmetric space (3.2) is the root lattice C n of the symplectic 
group Sp{2n) (see Table 1), with multiplicities mg i ±g j = 2, m,2e i = 1 for i,j = 1, . . . , n, i < 
j [25, 27]. The coset space measure (2.25) for the case (3 = 2 and an even number of fermion 
flavours is thereby given as 

n 

du 2 (U) = [dV] 11 d n sin 2r; J] sin 2 ^ - Tj ) sin 2 (r; + Tj ) (3.5) 

i=l i<j 

with ri G [0, |]. By using the trigonometric identity 

2 sin(rj — tj) sin(rj + rj) = cos 2r^ — cos 2rj (3.6) 
and defining Aj = cos2rj G [—1, 1], we may write (3.5) in the simpler form 



d^ 2 (C/) = ffl ndA»A(A 1 ,...,A„) 2 . (3.7) 



i=l 



For any X G u(n) represented as a diagonal n x n matrix with entries x« G E, we have introduced 
its Vandermonde determinant which is the Weyl determinant (2.15,2.18) for the unitary group 
G = U (n) given by 

&U(n)i x ) = 11 {xi ~ Xj) = A(xi,. . . ,x n ) . (3.8) 

i<j 

In the polar decomposition (2.24), we define the n x n diagonal matrix 

-C "J 

and write the radial coordinates as 

i? = exp ( i a 1 (g> p) (3.10) 

which, since {a 1 ,^} = 2<5 U , has the requisite property that it anticommutes with the parity 
matrix in (3.1) as 

RT = FR- 1 . (3.11) 
We decompose the angular degrees of freedom V G U(n) x U(n) as 

^=r + vl) (3 - i2) 



ii 



with V± G U(n), so that they commute with the parity matrix as 

VF = TV . 



(3.13) 



On substituting (2.24) into the finite volume partition function (2.6), and by using the identities 
(3.10)-(3.13) along with cos^ 1 ® 77) = I2 ® cost; for any 77 G u(n), we may thereby write the 
action of the non-linear cr-model in this case as 

Re Tr (MUFU^ = \ Tr [M F V (R 2 + iT 2 ) V' 1 ] 
= Tr [M F V cos (2a 1 ® p) V' 1 } 



Tr 



'Ai 



V, 



Tr 



Ml 



V+ 1 



V 



-1 



An , 



(3.14) 



The resulting integrals over the two sets of angular coordinates V± G U (n) are identical and 
can be evaluated by applying the Harish-Chandra formula (2.26) to the unitary group G = U(n) 
(In this case (2.26) is often refered to as the Itzykson-Zuber formula [37] in the physics literature). 
The Weyl group in this case is the symmetric group Wut n ) = S n and it acts on X by permuting 
its eigenvalues, {X w )i = (w X w~ 1 )i = x^^, where w G S n is a permutation and sgn(ti;) is its 
parity. In this way the finite volume partition function (2.6) acquires the form 



zf n) (/Ji,...,/j n ) 



A( W ,...,/^^ 



n 

sgn(^+tZ>_) J|' / dAi e A ^^ + W+^-W) (3.15) 
i=i •/, 



involving the unfolded masses /jj, i = 1, . . . ,n. Since the integration measure and domain in 
(3.15) are invariant under permutations of the Aj's, we may reduce the double sum over the 
Weyl group to a single sum over the relative permutation w = w + wZ 1 - The Aj integrals are 
elementary and expressed in terms of the function 

1 



K(a;) 



sinh x 1 



x 



-1 



dA e 



x A 



(3.16) 



which for purely imaginary argument defines the spectral sine-kernel of the unitary ensemble [32] . 
The sum over w G S n defines a determinant, and in this way one may finally arrive at the simple 
n x n determinant formula 



Zf n) (/ii,...,/z n ) 



det 



K(//j + fij) 



i,j=l,...,n 



A(m, . . . ,/i n ) 2 



(3.17) 



This is precisely the form of the finite volume partition function that was obtained in [19] from 
considerations based on partially-quenched effective field theories. 
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N f = 2n + l 



The restricted root lattice of the symmetric space (3.3) is the non-reduced irreducible root lattice 
BC n given by the union of the root systems B n and C n of the groups SO(2n + 1) and Sp(2n) 
(see Table 1), with multiplicities mg i ±g j = 2, ni2e i = 1 as above, and mg i = 2 [25, 27]. In this 
case the measure (2.25) thus reads 



du2{U) = [dV] JJdrj sin2rjsin 2 rj sin 2 (rj — rj) sin 2 (rj + rj 



i=i 



i<j 



[dV] 



2 n 



H dA * A(A 1; ...,A n ) 2 , 



(3.18) 



i=l 



where we have used the double angle identity sin 2 rj = \ (1 — cos 2rj). With the definition (3.9), 
the radial coordinate matrix 



R = 



exp ( i a 1 ® p) 



(3.19) 



satisfies the requisite condition (3.11). The angular matrices V G U(n+\) x U(n) are decomposed 
as in (3.12) with V + G U(n + 1) and F_ G U(n), and the non-linear c-model action is now 
modified to 



Re Tr 



Tr 



+ Tr 



Mo 



Mi 



Mn / 



/l 



V 



+ 



Ml 



\ 



An/ 



Mn ■ 



(3.20) 



The Harish-Chandra formula applied to the integrations over V+ G U(n + 1) and V- G U (n) 
then yields 



Z. 



(2n+l) 



(Ho, Hi, . . . ,// n ) 



jl (/*» — A*o) 1 

X? X] y sgn(u> + ) sgn(u>_) e MtS + (0) 

A( m ,..., Mn ) ^ 

x[] / dA, e^^+W+^-W 5 + (-1)^ ■ 



(3.21) 

Note how the extra factors of (1 — Aj) in (3.18) have nicely cancelled against the (n + 1) x (n + 1) 
Vandermonde determinant A(l, Ai, . . . , A n ) coming from the V+ integral. In the sum over u>_ G 
S n , we may regard S n as the subgroup of the permutation group S n +i acting only on the indices 
i = 1, . . . , n in {0,1, ... ,n}. Then as before we can truncate the double Weyl group sum in 
(3.21) to a single sum over w = w + wZ 1 G S n+ i by extending w_ to all of S n+ i through the 
definition w_ (0) = 0. After performing the Aj integrals, the finite volume partition function is 
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thus given by the simple (n + 1) x (n + 1) determinant formula 



fl (JH ~ W)) 1 . 

„(2n+l)/ \ i=l / 

A( W ,...,/i n ) 2 \ 

+ (-l) n+JV P det 



det 



g MO e A* j 

_K(^ + Mi) K(/Xi + /ij)_ 

e -Mo g — 

K(/i + /Xj) K(//j + //j) 



i,j=l,...,n 



i,j=l,...,n. 



(3.22) 



This is a new simplified expression for the finite volume partition function of QCD3 with an odd 
number of fundamental fermions. It can be readily checked to coincide with previous expressions 
(given by (2n + 1) x (2n + 1) determinant formulas [15, 20]) for low numbers of flavours. 

3.2 Correlation Functions 



For (3 = 2, the fc-level correlation functions can be expressed as a ratio of partition functions 
with Nf and Nf + 2k flavours as in (2.8). By using (3.17) one finds the scaled correlators for 
an even number of quark flavours given by 



det 



(2ra)/> > 

P2 Isi) ■ ■ ■ >sfc;Mi> • 



K(//j + /Xj) K(fj,i + i(i> 
K(^--iCi) K(i0-i0' 



i,i'=i,. 



det 



(3.23) 



where the matching condition (2.9) in this case fixes Cf } = (^) fe . This again agrees with the 
field theoretic results of [19]. In the form (3.17), the expression (2.8), originally derived in the 
context of random matrix theory [15], is thereby naturally tied to the field theory derivation of 
the spectral functions ■ 

For an odd number of flavours one finds 



(2n+l),^ A 

det 


in) 

e M e^i e' C i' 
K(/io+w) K(/if+/Xj) K(^+iC;/) 
.K( W) -iC() Kfo-iG) K(iCi-iCiO. 


+ (_i) fc +" +Jv p J w^-w L 

i,j=l,...,n I Ci^-Ci J 
«,/'=!,. ..,fc 


(2n) k det 


c MO c M J 


+ (-!)"+"* • 

i,j=l,...,n 


M0-+-MO 
Pi^-Pi 


\ 



(3.24) 

This is a new general expression for the spectral fc-point functions in the case of an odd number of 
quarks at (3 = 2. It coincides with the known expressions from random matrix theory [16] and the 
replica method approach to QCD3 [21]. In particular, the formula (3.24) shows explicitly that the 
correlators are not even functions of the scaled Dirac operator eigenvalues Q for M / O271+1 [21]. 
This expression also covers the case of a single flavour Nf = 1 (n = 0), which can be obtained 
by omitting the corresponding n x n blocks in the determinants above. 
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3.3 Sum Rules 



We may work out the massive spectral sum rules (2.13) by taking the equal mass limits in (3.17) 
and (3.22). This yields indeterminate forms which can be computed by regulating the mass 
matrices in any way that removes the n-fold and (n + l)-fold eigenvalue degeneracies, and then 
taking the degenerate limits using l'Hopital's rule. However, this is rather cumbersome to do in 
practise and it is much simpler to go back and work directly with the radial coset representation 
of the finite volume partition functions. 



Nf = 2n 



From (3.14) we see that in the case of equal quark masses \n = fi Vi = 1, . . . ,n, the u-model 
action is independent of the angular degrees of freedom V± £ U(n) and the low-energy dynamics 
is mediated entirely by the n radial Goldstone boson degrees of freedom. The partition function 
(2.6) is then simply given by 



n „ 

Zf n) ( M , • • • ,/i) = [J / dA, e 2 ^* A(A 1; . . . , A n ) 

i=l _-y 



(3.25) 



We now use the identity [32] 



A(Ai,...,A„) = det 



Ar 1 



i,j=l,...,n 



(3.26) 



and expand the resulting two determinants in (3.25) into sums over permutations w± £ S n . As in 
Section 3.1, the partition function depends only on the relative permutation w = w + wZ 1 £ S n , 
and the Aj integrals are now expressed in terms of derivatives of the spectral kernel defined by 
(3.16) as 



i 

K (m) ( ) = d^K(x) = 1 f dxx m eXX 

y ' dx m 2 J 

-1 



(3.27) 



In this way we arrive at the n x n determinant formula 

zf n) (/i,..., / u) = det[K^- 2 )(2 /U )' 



i,j=l,...,n 



(3.28) 



We may present the partition function (3.28) in a somewhat more explicit form by expressing 
the generic derivatives (3.27) in terms of the sine-kernel K(x) itself. Starting from the integral 
representation (3.27) we write 



J dX e xX I 1 - m J 



ds s 



m—l 



-1 



ds s m - 1 (e xs - e~ x ) 



(3.29) 



15 



This yields a simple linear inhomogeneous recursion relation for the sequence of functions 
K( m )(x), m > in (3.27) given by 



KM(x) + ^ K^" 1 ^) = K(x) + ^ (l - (-1) 



"' < < — KS, 



(3.30) 



It is straightforward to iterate the recurrence (3.30) in m and express the solution in terms of 
K(x) and other elementary functions of x, and in this way the equal mass partition function 
(3.28) can be written as 



Z^ n \(i, ...,/x) =det 



K 



i+j-2 



i,j=l,...,n 



(3.31) 



where 



K m (x) 



2 ^ (m- k)\ V x) 



k+l 



1 x - ml 

[(l - (-l) m+fc ) (smhx - coshx) - 2sinh: 



(3.32) 



From (3.31) one may now proceed to derive the massive spectral sum rules (2.13). For the 
first few even numbers of quark flavours we find 



C>° I 2 



c>° / 2 



AT /= 2 



iVf=4 



N f =6 



2/t coth 2/i— 1 
47? 



1 2sinh 2 2/j-/t sinh4/j-4/j 2 
4/77 4/t 2 -sinh 2 2/i 



x 2/ t (l6//+60/r 2 +3) cosh2/t-28/t 2 (4// 2 +3) sinh2/t-6/t cosh 3 2/t+9 sinh 3 2/x 
12/t 2 4// 2 (4/i 2 +3) sinh 2/t- 16/t 3 cosh 2/t-sinh 3 2/i 



These all agree with (2.14) in the limit /j, 
obtained in [14]. 



(3.33) 

0. The first two sum rules in (3.33) were also 



N f = 2n + 1 

From (3.18) and (3.20) it follows that in the equal mass limit /xo = Hi = H Vi = 1, . . . ,n the 
finite volume partition function (2.6) is given by 

n \ 

Zf n+l \^...^)= e" J] jdX^l-X,) e 2 ^ A » A(A 1 ,...,A n ) 2 + (— 1)^ {/x — > —fi} . 



(3.34) 
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Proceeding exactly as above, and using the fact that the functions (3.32) obey the reflection 
property K m (—x) = (— l) m K m (x), we can write (3.34) explicitly as the n x n determinant 
formula 



zf n+1) Gu, . . . = e" det[K i+j _ 2 (2fi) - K i+J _ 1 (2^) 



i,j=l,...,n 



+ (-l) JV P e"" det K i+j _ 2 (2fi) + Ki+j^fi) 



i,j=l,...,n 



(3.35) 



with the determinants omitted for the n = case as before. For an even number Np of Dirac 
operator eigenvalues, from (3.4) and (3.35) we obtain the massive spectral sum rules 



§< 2 +"7 2 



Nt=l 



N f =3 



tanh /j, 
2/j, 



1 2/j, 2 coth/x — + 3cosh2^) — 2sinh2^ 
6/i 2 2^i + sinh2/i 



(3.36) 



which again agree with (2.14) in the limit /i — > 0. 



4 Pseudo-Real Fermions 

The next simplest case is the parity-invariant Dirac operator i Jp acting on iV c = 2 colours of Dirac 
quarks in the fundamental representation of the SU(N C ) gauge group. In contrast to the previous 
case, the Dirac operator now possesses a special anti-unitary symmetry, [€, \If)\ = 0, generated 
by an operator (£ with (£ 2 = l [9, 10]. The operator £ can be expressed as a combination of 
the usual charge conjugation C = \a 2 with C (a^ 1 )* C^ 1 = cr M for n = 1,2,3, an infinitesimal 
rotation by the matrix id 2 in colour space, and complex conjugation. In a basis of spinors IV'fc) 
wherein <t\ipk) = l^fc), the matrix elements (i/jk\ i-P l^z) are all real [10]. This basis is obtained 
by augmenting the flavour components of the fermion fields by the two-colour group, thereby 
enhancing the original U{Nf) flavour symmetry to a higher pseudo-real global flavour symmetry 
group Sp(2Nf). In matrix model language this corresponds to the orthogonal ensemble (3=1. 

The group USp(2n) acts on the vector space C 2n preserving its canonical inner product as 
well as a symplectic structure $ n : M? n — ► M 2n , i.e. a real antisymmetric 2n x 2n matrix with 
J 2 = — l 2ra - In a suitable basis of C 2n , it consists of unitary matrices V satisfying 

VS n V T = J n , J n = ia 2 0lL n . (4.1) 
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With respect to the same basis, the mass and parity matrices are given by 



M 



a 3 <g> i a 2 



( Mo i o- 2 
V 



Ml 



a 3 <8) i a 2 



) 



N f = 2n, 



N f = 2n+1 



a 3 <S> In 



1 a 



a" 



N f = 2n, 
N f = 2n + 1 



(4.2) 



The antisymmetric forms of these matrices follow from the Grassmann nature of the quark 
spinor fields ipk described above. Again the mass matrices break the U Sp{4n) and U Sp(4n + 2) 
flavour symmetries, and the corresponding Goldstone manifolds are 



& (2n) = USp(4n) / USp(2n) x USp{2n) , 
<?i (2n + 1) = ?75p(4n + 2) / ?75p(2n + 2) x USp(2n) 



(4.3) 
(4.4) 



of real dimensions dim(?i(2n) = 4n 2 and dim^i(2n + 1) = 4n (n + 1). As expected, the number 
of Goldstone bosons here is exactly twice that of the previous section. 



4.1 Partition Functions 



N f = 2n 



The restricted root lattice of the symmetric space (4.3) is again the root space C n of Sp(2n), 
but now with multiplicities mg i ±g j = 4, m2g i = 3 [25, 27]. The corresponding radial integration 
measure (2.25) on Qi(2n) is thus 



du^U) = [dV] Yldn sin 3 2rj sin 4 (ri - Tj ) sin 4 (ri + rj ) 



i=i 



= JSiyndA^l-Af) A ( Ai,...,A n) 4 , 



2 n(2n-l) J 



(4.5) 



i=l 



where we have again defined Aj = cos2rj £ [—1, 1]. With p defined as in (3.9), we parametrize 
the radial element R as 



i? = exp (ia 1 <g> £) 

where the matrix 

i = a 2 ®p 

commutes with the symplectic structure as 



(4.6) 



(4.7) 



(4.8) 
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It follows that .R T = R^ 1 = w, and the radial matrix (4.6) anticommutes with the parity 
matrix T as in (3.11). For the angular degrees of freedom V G USp(2n) x USp(2n) we again 
take the decomposition (3.12) with V± G USp(2n), so that from (4.1) we have 



yry 



r . 



(4.9) 



Using these relations, along with the identity cos(<7 2 <8> rj) = I2 <8> cos?? for any 77 G u(n), we 
can now write the action of the non-linear er-model (2.6) similarly to (3.14) as 



Re Tr [MUYU T ^j = \ Tr [TMV (R 2 + R~ 2 ) V' 1 } 
= Tr [r M V (1 4 ® cos 2/?) V^ 1 ] 



= Tr 



mi 



+ Tr 



Mi 



lb 





(4.10) 



The finite volume partition function (2.6) thereby becomes 

1 



lb n 

Z[ 2n \ t i 1 ,...,» n ) = l[ dX, (1 - Af ) A(A 1 ,...,A n ) 4 0S n) (A 1 ,...,A n ;/x 1 ,...,/x n ) 2 (4.11) 
i=i ^ 



where 



^1 ^(Ai, . . . , A n ; fii, . . . , fi n ) 



J [dV] exp Tr 



Sp(2n) 



Ml 




V 1 2 




(4.12) 



It is not at all clear at this stage how to proceed with the calculation. The Lie algebra sp(2n) 
consists of 2n x 2n matrices X obeying the equation X J n + J„ X T = 0. The diagonal source ma- 
trices in (4.12) are not elements of this Lie algebra, nor are the original mass and parity matrices 
in (4.2) elements of sp(4n). Thus the Harish-Chandra formula (2.26) cannot be immediately 
applied to the angular integral (4.12). For the same reason this integral cannot be evaluated by 
using character expansion techniques for the symplectic group Sp(2n) [38]. It would be inter- 
esting to determine the corresponding analytical continuation of the Harish-Chandra formula to 
this case, in analogy with that done in the unitary case G = U (n) [39] which is based on the fact 
that the representation theories of the Lie groups U(n) and GL(n,C) are essentially the same. 
Note that the integral n^(A; fi) can be computed by means of the Harish-Chandra formula for 
SU{2) Sp{2). 
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In [18] the finite volume partition function in this case is evaluated by applying a formalism 
based on skew-orthogonal polynomials in the orthogonal ensemble of random matrix theory, 
resulting in the Pfaffian expression 



a,b=l,...,2n 



PfafF K' (fl a - Aft) 

. . . , Mn ) = i Wx,„.,z» (413) 

n Mi A(/xf,...,/x2) 2 
1=1 

where M2i-i = — Mi) M2i = Mi f° r * = 1, ■ ■ ■ Using the permutation expansion of the Pfaf- 
fian, the integral representation (3.27) and by comparing (4.13) with (4.11), it is possible to 
write a conjectural formula for the integral (4.12). There are several qualitative features of the 
expression (4.13) that make its derivation directly in field theory appear feasible. First of all, 
the denominator of (4.13) can be expressed in terms of the square of the Weyl determinant 
Ag p ( 2n )(X) of an element X G sp(2n) with eigenvalues Mi> . . . , [i n (see (2.15) and Table 1), as 
appears in the expression (2.26). Furthermore, the Weyl group W Sp ( 2n ) is composed of permu- 
tations of eigenvalues x% i > x n ^ , i — 1, . . . ti, 7r G S r ij along with, reflections x% i > x%. After the 
appropriate radial integration, the sum over Ws p (2n) can thereby produce the required Pfaffian 
in the numerator of the expression (4.13), which in the random matrix theory calculation arises 
naturally from quaternion determinants generated in the skew-orthogonal polynomial formal- 
ism. In the next subsection shall explicitly work out the <r-model partition functions in the case 
of equal quark masses and thus provide further evidence that the direct field theory calculation 
indeed does reproduce this Pfaffian expression. 



N f = 2n + 1 

The restricted root lattice of the symmetric space (4.4) is again the root space BC n , with 
multiplicities mg.±g. = 4, m 2 g x = 3 as above, and in addition mg. = 4 [25, 27]. The coset space 
measure (2.25) in this case is thus given by 

n 

dvi(U) = [dV] JJdrj sin 3 2r^ sin 4 rj sin 4 (r^ — rj ) sin 4 (r« + rj ) 

i=l i<j 

= ^Sri) fl dA * C 1 "*?) C 1 "**) 2 ACAi,...,^) 4 . (4.14) 

i=i 

With the definitions (3.9) and (4.7), the radial element 



<7 2 



R= { exp(i^ {) ) <415) 
satisfies the required condition (3.11) and also R T = R^ 1 = R' . The angular matrices V G 
USp{2n + 2) x USp{2n) are decomposed as in (3.12) with V+ G USp(2n + 2) and V- G USp(2n), 
so that (4.9) is obeyed. The cr-model action can be worked out analogously to (3.20) and (4.10), 
giving 

l 

\4 



11 r. 

z} 2n+1) ( Mo ,Mi,---,Mn) = II MMl-A 2 ) (1-A,) 2 A(Ai,...,A r 

i=l i 



r 



x Ai, . . . , A n ; Mo, Mi) ■ ■ ■ ) Mm) ^^(Ai, • • • , A n ; fx±, . . . , fj, n ) 

+ (-1)^ {%ZZX} . (4.16) 
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For n = this expression is understood as being equal to £l\ (l;/xo) + (—1) ^ Q\ (1; — A*o)- 
In [18] an expression for the finite volume partition function in this case is derived from random 
matrix theory in the case of a vanishing unpaired quark mass = (so that the matrix model 
measure is positive and the partition function vanishes when Np is odd). In the next subsection 
we derive a formula for the equal mass limit of (4.16) without this restriction. 



4.2 Sum Rules 



N f = In 

As in the previous section, to work out the massive spectral sum rules (2.13) in the even-flavoured 
case at (3 = 1 it is convenient to work directly with the partition function (4.11) for equal quark 
masses [i{ = fi Vi = 1, . . . ,n, which is determined entirely by the n radial Goldstone boson 
degrees of freedom and is simply given by 

n \ 

zfV...,/i) = ;Q /dA 4 (l-A 2 ) e 4 " A * A(A 1; ...,A n ) 4 . (4.17) 
i=i \ 

Because of the original fermion determinant required in (2.12) and the doubling of the flavour 
symmetry in this case, we can anticipate the appearence of a Pfaffian expression for the finite 
volume partition function, in contrast to the determinant formula of the previous section. For 
this, we write the power of the Vandermonde determinant in (4.17) as the determinant of a 
2n x 2n matrix given by [32, 40] 

A(A 1 ,...,A n ) 4 = det[ A | 1 , 2 j i=Wl . (4.18) 

LU Jj=l,...,2n 

Expanding the determinant into a sum over permutations w £ S^n, and using the permutation 
symmetry of the integration measure and domain in (4.17) to appropriately rearrange rows, we 
may bring the equal mass partition function into the form 



zfv. ..,/,) = y, s ^ n ^ n / dA * t 1 -^ 2 ) 

- =s 2 „ i=l _i 

(w(2i) - l) . (4.19) 



X 



Let S 1 ^ denote the subgroup of S^n consisting of "increasing" permutations w, i.e. those obeying 
w(2i) > w(2i — 1) Vi = 1, . . . , n. Then the expression (4.19) can be written as 



Zi 2n \»,...,») = Yl ft /dA, (l-Af) 

wPS> 1=1 i 



tses>, i=i _i 

x (w(2i) - w(2i - 1)) A f (2i)+ " (2l - 1} - 3 . (4.20) 

The expression (4.20) now has the standard form of a Pfaffian [32] and we may thereby write 
it as 

Zf n) = Pfaff A 2n (4fi) , (4.21) 
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3 giA 



where A2 n {x) is the In x 2n antisymmetric matrix with elements 

1 

A 2n (x) ij = (i - j) J dX (1 - A 2 ) A i+ ' 
-l 

= 2(i - j) (K^'- 3 )(x) - K^'" 1 ^)) 
These matrix elements can be alternatively expressed as 



A 2n {x)ij = 2(i - j) 



d 
dx 



i+j-3 



ds s I dX e xX = {i-j) 



dx 



K'(x) 



(4.22) 



(4.23) 



The substitution of (4.23) into (4.21) can be straightforwardly shown to coincide with the equal 
mass limit of the finite volume partition function (4.13) obtained from the orthogonal ensemble 
of random matrix theory. 

Proceeding as in Section 3.3 to rewrite derivatives of the spectral sine-kernel in terms of the 
functions (3.32), the equal mass partition function in the case of an even number of pseudo-real 
quarks acquires the Pfaffian form 



(4.24) 




It can be computed explicitly by using the Laplace expansion of the Pfaffian [32] 



2n 

Pfaff A = J2 

3=1 



(A)ij Pfaff A 



(4.25) 



for any fixed i = 1, . . . , 2n, where the Pfaffian minor Ajiji is obtained from the antisymmetric 
2n x 2n matrix A by deleting both its i-th and j-th rows and columns. The massive spectral 
sum rules can thus be extracted from (3.32) and (4.24) by the substitution Nj — > 2Nf in (2.13) 
to account for the symplectic doubling of fermion degrees of freedom in this case, and for the 
first few even numbers of quark flavours we find 



i 



N t =2 



N f =A 



Aii coth 4/i- 1 16/? 



I 36/i(4/i 2 +5) sinh8At-9(64/t 2 +5)sinh 2 4/^-16/i 2 (256/t 4 +216/t 2 +27) 
16m 2 9(8Ai 2 +l)sinh 2 4/t-36/i sinh8/i+128/i 6 (16/i 2 +9) 



(4.26) 



These results agree with (2.14) in the limit u — > of massless quarks. 
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N f = 2n + l 



The equal mass limit Hi = H Vz = 0, 1, . . . , n of the partition function (4.16) is given by 



n ,. 

Zf n+l \n,...^) = e^I] jd\ t (l-\i) (1-A,) 2 e 4 ^* A(X 1 ,...,X n y 



+ (-1)^ 



(4.27) 



The calculation is identical to that of the even-flavoured case above. The only difference is that 
the antisymmetric matrix (4.22) now contains an additional measure factor (1 — A) 2 , which when 
expanded leads to the equal mass finite volume partition function in the case of an odd number 
of symplectic fermions in the Pfaffian form 



Z< 2n+1 V-^)=e 2 ^ Pfaff 



(i-j) K l+j _ 3 (4 At )-2K l+j _ 2 (4/ i )+2K l+J (4 M )-K l+j+ i(4 M ) 



+ (-l) N P c~ 2 » Pfaff 



i,j = l, . . . ,2n 

(i-j) (K l+J _ 3 (4 M )+2K l+J _ 2 (4 At )-2K l+j (4 M )-K, + , +1 (4^) 



i,j=l,...,2n 



(4.28) 

This expression also covers the case Nf = 1 of a single fermion of mass /i, obtained by dropping 
the Pfafnans in (4.28), which for an even number Np of Dirac operator eigenvalues leads to the 
massive spectral sum rule 



(4.29) 




Sum rules for higher numbers of flavours are also straightforward to work out, but they 
become increasingly complicated and will not be displayed here. In each case one finds that the 
massless spectral sum rule (2.14) for (3 = 1 is indeed obeyed. The equations (2.10), (2.11), (4.24) 
and (4.28) suffice to express to smallest eigenvalue distribution of the QCD3 Dirac operator for 
SU{2) colour group in a closed form for an arbitrary number of fermion flavours in the case of 
equal quark masses. 



5 Adjoint Fermions 

Finally, we consider the case of quark spinors which transform in the adjoint representation of 
the SU{N C ) colour group for N c > 2, which is real. As in the previous section, the adjoint 
representation of the Euclidean Dirac operator i ]p uniquely defines an anti-unitary operator £0 
with [£o,i-^] = 0, but now 

(£ ) 2 = -1 [9, 11]. The operator £0 is given as a combination of 
charge conjugation C = iu 2 and complex conjugation. In this case the Dirac operator can be 
expressed in terms of real quaternions [11], and the condition (tolV'fc) = iV'fc) defines a basis \ipk) 
of Majorana spinors of the second kind. The generic unitary U(Nf) flavour symmetry thereby 
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truncates to its real orthogonal subgroup O(Nf). In matrix model language this corresponds to 
the symplectic ensemble (3 = 4. The mass and parity matrices are given exactly as in (3.1), and 
so the Goldstone manifolds are 

G 4 (2n) = 0(2n) / 0{n) x O(n) , (5.1) 
G A {2n + l) = 0(2n + l)/0(n + 1) x O(n) (5.2) 

of real dimensions dim^4(2n) = n 2 and dimQ A (2n + 1) = n (n + 1). As expected, the number 
of Goldstone bosons here is exactly half that of Section 3. 



5.1 Partition Functions 



N* = 2 



The "degenerate" case Nf = 2 corresponding to two-flavour component adjoint fermions is very 
special because it represents flavour symmetry breaking from 0(2) to a finite subgroup Z 2 x Z 2 . 
The coset space (5.1) in this case is the orbifold 



&(2) = 0(2) /Z 2 x Z 2 . 



(5.3) 



The quotient by one of the Z 2 factors in (5.3) can be used to eliminate one of the connected 
components of the 0(2) group, thereby realizing the orbifold explicitly as a one-dimensional 
manifold diffeomorphic to an open interval 



g A {2) = 50(2)/Z 2 * (0,tt) 



(5.4) 



After flavour symmetry breaking, there is thus a single massless degree of freedom described in 
the mesoscopic scaling regime by a real scalar field 6 G (0,7r). 

To evaluate the corresponding finite volume partition function (2.6), we parametrize the 
elements U G 50(2) appearing in (5.4) as 



U 



COS I 

sin( 



- sin( 
cos 6 



and work out the effective cr-model action to be 

Re Tr (m UT U t ) = 2/i cos 26 



(5.5) 



(5.6) 



With the normalized Haar measure [dU] = d6/Tr on 50(2), the partition function (2.6) in the 
two flavour case is thereby found to be 



zi 2 V) = I„(2p) 



(5.7) 



where generally 



l{x) = J 



, x cos <j>+i m < 



(5.8) 



is the modified Bessel function of the first kind of integer order m. 
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N f = 2n 



Let us now turn to the case of generic even integer values Nf = 2n with n > 1. The restricted 
root lattice of the symmetric space (5.1) is the root space D n of the special orthogonal group 
SO(2n) (see Table 1), with root multiplicities mg i ±g j = 1 for i,j = l,...,n,i< j [25, 27]. The 
corresponding radial integration measure (2.25) on Qi(2n) is thus 



dv±{U) = [dV] JJdrj Y\ | s i n ( r i _ r j) sin(r; +Tj) 

i<j 



i=l 



(5.9) 



Similarly to the n = 1 case, the quotient in (5.1) can be used to eliminate the connected 
components of the orthogonal groups, at the expense of presenting the coset space as the orbifold 



Q A {2n) = SO(2n)/SO(n) x SO{n) x Z 2 . 



(5.10) 



In this case, it is more convenient to parametrize the radial directions of this orbifold in terms 
of the angle variables 8{ = 2ri G [0,7r], and using (3.6) we may write (5.9) as 



du 4 (U)= } dV [ u UdOi A(cos0i,...,cos0 n ) 
22 n ( n + 1 ) ■ LJ - 



(5.11) 



To ensure positivity of this measure when the absolute value signs are dropped, we will restrict 
the integration region for the di's to the domain 



{(01,... A) | O<0!<0 2 <...<0 n <7r} 



(5.12) 



This can always be achieved without complications for the integration of symmetric functions. 

The particular details of the radial decomposition here are identical to those of Section 3.1, 
and we thereby arrive at the non-linear a-model action 



Re Tr {mUTU^^ = Tr 



mi 



' COS 01 



v_ 



+ Tr 



Ml 



+ 



V 



V 



-1 



cos n 



' COS 01 



V 



-1 



cos 6 n 



(5.13) 



where now V± £ 0(n). The finite volume partition function is thus given by 

/n 
'[[dOi A(cos6>i, . . . ,cos0 n ) f4 n) (cos6>i, . . . ,cos0 n ;/xi, . . . ,ji n ) 2 (5.14) 



i=l 



where 



^4 ^ ( x l , • • • j x n i Ml 5 • • • ) Mn) 



J [dV] exp Tr 



XI 



V 



(5.15) 



SO(n) 
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As in the previous section, the integral (5.15) cannot be evaluated using the Harish-Chandra 
formula (2.26) because the Lie algebra so(n) consists of antisymmetric nxn matrices. The first 
two members of this sequence of angular integrals are easily evaluated explicitly to be 

n£\x 1 ,x 2 ; l i 1 ,H2) = 2e5("^X^>/ (i( ft - M2 )(x 1 -i 2 )) . (5.16) 

For n = 3 the integral (5.15) can be computed by means of the Harish-Chandra formula for 
SU{2) 50(3). 

In [18] the finite volume partition function in this instance is computed using a skew- 
orthogonal polynomial formalism in the symplectic ensemble of random matrix theory. For 
n = 2£ with I > 1 it is necessary to assume that the mass eigenvalues \i{ appear in degenerate 
pairs, i.e. fj,2a-i = l^2a for each a = 1, . . . ,£, in order to ensure positivity of the matrix model 
measure. In this case a PfafRan expression for the partition function is obtained in terms of 
integrals of the spectral sine kernel (3.16) (rather than derivatives as in the previous section), 
or equivalently in terms of exponential integral functions Ei(x). In general, this restriction does 
not aid the field theory calculation, which should thereby capture the case of generic mass eigen- 
value configurations. In the next subsection we will evaluate (5.14) in the equal mass limit and 
find a striking contrast between the field theory and random matrix theory expressions. As in 
(5.7), the field theory calculation will naturally involve modified Bessel functions of the first kind 
(along with sine kernel derivatives) , which is reminescent of results for QCD in four spacetime 
dimensions [2, 30]. This may be a consequence of the known relations between the QCD3 and 
QCD4 partition functions in the e-regime [20, 22]. Comparison of our expressions with those 
of [18] imply that some non-trivial simplifying identities among these Bessel functions ought to 
hold, but we have not been able to find any. 

For n = 21 + 1, in addition to pairwise degeneracy, it is necessary to assume in the ran- 
dom matrix theory approach that the single unpaired quark mass vanishes, again for positivity 
reasons. In the next subsection we derive the equal mass limit of (5.14) without this restric- 
tion. From the field theory perspective, we shall also encounter a distinction between the cases 
Nf = A£ and Nf = A£ + 2, and this distinction is very natural. From Table 1 we see that the 
Weyl group sum for the group SO{2£ + 1) is the same as that of the symplectic group Sp(2£). 
In contrast, the Weyl group of SO(2£) contains the usual permutation part Si and the factor 
(Z2)^ _1 is the normal subgroup given by the possible assignments of an even number of sign 
flips to the eigenvalues. These facts, along with observations similar to those spelled out in 
the previous section, support the precise matching between the field theory and random matrix 
theory results in these instances. 

N f = 2n + 1 

The restricted root lattice of the Goldstone manifold (5.2) is the root lattice B n of the special 
orthogonal group SO(2n + 1) (see Table 1), with multiplicities mg i ±g j = mg i = 1 for i,j = 
1, . . . ,n, i < j [25, 27]. The corresponding coset space measure (2.25) is thus 
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dv^U) = [dV] Y\_dr-i sinrj sin(rj — rj) sin(rj + Tj) 



i=i 

n 



[dV] l\d Ki I A («?,..., 4)| , 



(5.17) 



i=i 



where we have defined Kj = cosrj G [0, 1] and used (3.6) along with the double angle identity. 
By introducing the domain 

K = {(«!,..., K„) | 0< Ki < K 2 < ... < K n < l} (5.18) 

and proceeding exactly as in Section 3.1, the finite volume partition function may then be written 



as 



/n 
Y\ dK, t A (kI, ■ ■ ■ , kI) ^i n) (2k? - 1, . . . , 2/4 - . . . 



/in) 



x ^ n+1) (1,2k? - 1, . . . ,2k* - l;/u ,/ii, . . . ,n n ) 
+ (-1)^ {^=£°} 



(5.19) 



with the same conventions for the n = case as in (4.16). The odd-flavoured cases are not 
covered by the random matrix theory analysis of [18]. In the next subsection we will compute 
(5.19) explicitly in the case of equal quark masses. We will find that, in contrast to the previous 
two classes, for adjoint fermions they involve fundamentally different functional structures than 
those of the even- flavoured cases. In this case the partition function is expressed in terms of a 
combination of error functions and Hermite polynomials, and appears to be more akin to the 
random matrix theory predictions [18, 32]. These functional forms are similiar to those which 
arise in partially-quenched effective field theories of QCD3 [19]. 

5.2 Sum Rules 

To compute the massive spectral sum rules (2.13), we proceed as previously to compute directly 
the finite volume partition functions (2.6) in the limit of equal quark masses for the various 
congruence classes modulo 4 of the flavour number Nf. 



N f = 2 

For Nf = 2, we substitute (5.7) into (2.13) to compute 




(5.20) 



Once again this agrees with (2.14) in the limit fi — > 0. 
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N f = U 



For Nf = 2n with n > 1 the partition function (5.14) for equal quark masses is given by 

/n 
Hd9i e 2 ^ cos9 ' A(cos0i,...,cos0 n ) . (5.21) 

Again, since the functional integration in the original QCD3 partition function is over real 
Majorana fermion fields, thereby producing a square-root fermion determinant in (2.12), we can 
anticipate a Pfaffian expression for (5.21). With this in mind, we will proceed to evaluate the 
finite volume partition function using the method of integration over alternate variables [32, 40], 
whose details again depend crucially on the even/odd parity of the flavour rank n. 

For n = 21 with I > 1 we can insert the determinant representation (3.26) to write 



£ ® 2a 

zfV...,/*) = J f[d6 2a e 2 " 008 *- I d0 2a _! e 2 " 008 *- 1 detfcos'- 1 ^]. j=h 



.,21 



(5.22) 

where we have defined #0 = 0. We formally insert the integrals into the determinant in (5.22) 
by using its multilinearity to write 



1 

Zf\^...^)= / l\de 2a e 2 ^ cos ^ det 



cos^ 1 Q. 



02b 



2 b 



J d6 e 2 » cosd cosine 



(5.23) 



j=l,...,2£ 



Since the integrand of (5.23) is a symmetric function of the integration variables, we may extend 
the lower limit of integration of each variable inside the determinant down to 0. This symmetry 
also allows us to extend the range of each of the remaining integration variables to [0, tt]. After 
rearranging rows in the determinant and expanding it as a sum over permutations, we arrive at 
the expression 



Zf\^...,ix) = £ sgn(«>) fl Jd6 2a e 

02a 

J d6 e 2 ^ cose cos^ 20-1 )- 1 9 



2^cos6»2a cog «j(2a)-l q 



(5.24) 



which may be written as a sum over the subgroup C S 2£ of increasing permutations given 
by 



£ 7T 02a 

n) = ^2 sgn(w) f[ J d6 2a e 2 ^ cose2a J d9 e 2 ^ cosd 
weS> a=i 



^ cos ^(2a)-l cog »(2a-l)-l q _ CQS »(2a-l)-l ^ CQS ^(2a)-l ^ _ 

(5.25) 
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As before, the expression (5.25) now has the standard form of a Pfaffian and we may write 
it as 

Zf ] {p,,...,n) = Pfaff B 2l {2n) , (5.26) 
where B2t(x) is the 21 x 2£ antisymmetric matrix with elements 

7T 7T 

B 2 e(x)ij = J d6 e xcosd cos^H J d8' sgn (#-#') e xcosd ' cos^V . (5.27) 



The integrals in (5.27) can be evaluated explicitly by using the Fourier series representation of 
the sign function 



sgn(0) = "E 2fc + 1 ( 5 - 28 ) 
fc=0 

for 4> S [— 7r, 7r], and defining the sequence of functions 

7T 

l 2fc+ i(x) = J d6 e ^ose+(2fe+i)ie (5>29) 
o 

for I; £ No. We may then express the matrix elements (5.27) as series in derivatives of the 
functions (5.29) given by 

= ^ £ • (5 - 30) 

k=0 

Let us now explicitly evaluate the derivatives of the functions (5.29). The real parts can be 
written in terms of the modified Bessel functions (5.8) as 

Re\^l(x)=nlg- + 1 )(x) . (5.31) 

For the imaginary parts, we use the expansion 

sin(2£: + l)# = Im (cos 9 + i sm6) 2k+1 



sin 



»±Qtl)tyrQ cos*-*- e (5 .32) 

along with the change of integration variable A = cos# € [—1,1] in (5.29) to evaluate them in 
terms of the sine kernel derivatives (3.27) as 

Im l£T>) = 2 £ (f+l) £(-!)• Q KV-^Hx) ■ (5-33) 

The sums over I and a in (5.33) can be interchanged, after which the sum over p = I — a for 
fixed a can be evaluated explicitly by using the combinatorial identity 



p=0 
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that is straightforwardly proven by induction. Substituting in the explicit functional forms 
(3.32) for the sine kernel derivatives, and plugging (5.31) and (5.33) into (5.30) then enables us 
to write the equal mass partition function (5.26) explicitly as 



z! 4 V,...,/i) = Pfaff 



E 

k=0 



2k + 1 



i,j=l,...,2t 



(5.35) 



where 



r , , 1 A {-lf2 k - 2a (k-2a + i)\ (k-a\ ( _1\ 6+1 

*h;iW - ^ 2^ {k _2a-b + i)\ a i 

x I ( 1 — (— l) b ) ( sinhx — coshx ) — 2sinhx 



[(l - (-If) (sii 



(5.36) 



From (5.35) one may now extract the massive spectral sum rules (2.13) by using the derivative 
identities 



dl m {x) m dl k .i(x) 

dx = I m+l(x) + — I m (x) , ^ = I k ;i+l{X) 

For instance, for four flavours of adjoint fermions one finds 



(5.37) 



E 

fc=0 



2k+l 



N f =4 



4^ E 

fc=0 



J 2fc+l( 2 M) (l fc;1 (2 M )-^±ll fc . (2 M ))-.f 2fc+2 (2 M )I fc . (2 M ) 



2k + l 



(5.38) 



A numerical evaluation of (5.38) in the limit fj, — > establishes precise agreement with (2.14) in 
this case. 



Nf = 41 + 2 

Let us now consider the case Nf = 2n with generic odd flavour rank n = 2£ + l,^>l. The 
calculation proceeds analogously to that of the even rank case, except that now one needs to keep 
careful track of the single unpaired integration. Starting from (5.21) we derive the determinant 
formula 
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J f[d02a e 2 » cosd2a J d6 2a -i e 2 ^ 0892 ^ 1 

a=l n 

+1 02a - 2 



u 2t+l 



621 



i[d9 2a 



, 2/U cos 6» 2a 



det 



a=l 



2^+1 



02b 



COS-' 1 6*26 



/ d6> e 2 ^ cose cos^fl 

26-2 

J dcf) e 2 ^ cos * cos*' -1 



6=1,.../ 

j=i,...,2e+i 

(5.39) 



Using symmetry of the integrand in (5.39) to extend the integration ranges as before, expansion 
of the determinant into a sum over permutations yields 



W&S21 + 1 0=1 Q 

e 2a 



2Mcos6» 2a cog ?t>(2a)-l g 



X 



y d g e 2 M cose cog ri,(2a-l)-l ^ y d( £ e 2 M cos0 ^(2^+1)-!^ 



sgn(t&) ]J j d8 2a e 2fMCOs92a Jd8e 2 » cosd 
^es> +1 0=1 

(cos^)- 1 ^ cos^-^-^-cos-^- 1 )- 1 ^ cos*^)- 1 ^ 

(f) (5.40) 



X 



x y d(p e 




2m cos ^ cog ri)(2^+l)-l 



which again has the standard form of a Pfaffian of dimension (2£ + 2) x (2^+2). The last integral 
in (5.40) can be expressed in terms of derivatives of the modified Bessel function Iq(2^). With 
the same functions (5.36), we thereby find the equal mass finite volume partition function in the 
explicit form 





[ ~ 4- 1 1 ) (2M)I fe;J _i(2 M )-/0; 1 1 )(2M)I fc ; l -l(2M) J 
2^ 2k+l J (.-W 
fc=0 




= Pfaff 






-4 j '" 1) (2m) 


i,j=i,...,2e+i 



(5.41) 

The (rather cumbersome) sum rules may now be extracted from (5.41) exactly as described 
above for the even rank case. 
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N f = 4£ + 1 

For odd numbers of flavours Nf = 2n + 1 the equal mass limit of the finite volume partition 
function (5.19) is given by 



i=i 



P 1 :, 



(5.42) 



The calculation of the integrals in (5.42) is identical to that of (5.21), with the obvious modifi- 
cation of the integration domain and the replacements cos 6 { — > 2k? everywhere in the previous 
formulas. For n = 2£, one may in this way bring (5.42) into the Pfafnan form 

Zf +1) ( M , ...,»)= e-^- 1 ^ Pfaff E 2£ (4fi) + (-1)"* e^" 1 )" Pfaff E u (-An) (5.43) 

where Eitix) is the It x 2£ antisymmetric matrix defined by the elements 



E 



21 



;k' 2 {.J2\j~ 1 



( x ) {j = JdK e XK2 {k 2 ) 1 1 J dn' sgn(K-K') e XK " (k' 2 )' 





(5.44) 



In this case it is more convenient to use the Fourier integral representation of the sign function 
given by 



t \ 1 f i sin qy 
sgn(y) = - dq — — 
vr y q + ie 

— oo 



(5.45) 



for |/6R, where the infinitesimal parameter e — > + regulates the ambiguity in sgn(y) at y = 0. 
Then the matrix elements (5.44) can be written in terms of the functions 



E q (x) = J 



dn e XK +iqK 



(5.46) 



which may be integrated explicitly in terms of error functions. 

In this we may write the finite volume partition function (5.43) explicitly as 











= e"( 4 


'-V* Pfaff 


4tt J Ut * q+ie 

— oo 


i,,j=l,...,2* 


+ (- 


-l) N P {n - 







where 



(5.47) 



(5.48) 
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The matrix elements in (5.47) can be made somewhat more explicit by using the derivative 
identity 



d fc erf(x) _ . . k _, 2 



-2 



dxk = (-!)*"' e"* (5-49) 

where H m {x) are the Hermite polynomials. As before, this expression includes the one- flavour 
case Nf = 1 (I = 0), which is obtained from (5.47) by omitting the Pfaffians. These relations 
can also be used to now extract the massive spectral sum rules in these instances. 



N f = U + 3 

The partition function (5.42) for n = 21 + 1 is evaluated identically as in (5.40) with the obvious 
change of integration domain and the substitutions cos 6 i — > 2k? as above. The single unpaired 
integration in this case can be expressed in terms of derivatives of the function 

E o(x) = j- erf (iy^) , (5.50) 

and we can thereby write the equal mass finite volume partition function as the (2^ + 2) x (2^ + 2) 
Pfaffian formula 









= e -(«+i) n Pfaff 


— oc 






-S£ _1) (4ri 


i,j=l,...,2l+l 


+ (-i)^ w 







(5.51) 

This expression encompasses the case of three flavours (£ = 0) by omitting the (2£+ 1) x (21+ 1) 
blocks in (5.51). The Goldstone manifold (5.1) in this case is diffeomorphic to a disk, £4(3) — 
SO(3)/SO(2) x Z 2 = S 2 /Z 2 = D 2 , and the partition function is equal to the function (5.50) 
evaluated at x = 4/x. This simple evaluation follows of course from the explicit forms of the 
angular integrals (5.16) in this case. The corresponding massive spectral sum rule (2.13) is given 
by 




(5.52) 



In addition to the sum rules, the expressions (2.8), (5.7), (5.35), (5.41), (5.47) and (5.51) 
suffice to determine the spectral /c-point correlation functions in closed form for an arbitrary 
number of flavours Nf in the limit where the collection of equal quark mass and (imaginary) Dirac 
operator eigenvalues is completely degenerate. In particular, they may be used to determine 
the correlators at the spectral origin for massless quarks, as well as the smallest eigenvalue 
distribution functions (2.10,2.11). 
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